function OpenCRTBP_SM_BaseFigure
%% Constants
G           = 6.674e-11 * ((1/1000)^3);
mass_Mars   = 6.39e23;
mass_Deimos = 1.8e15;
mass_Phobos = 10.8e15;
GM_Mars     = G*mass_Mars;
GM_Deimos   = G*mass_Deimos;
GM_Phobos   = G*mass_Phobos;
GM_sun      = 1.327124400*10^(11);     
%[km^3/s^2] Gravitational Parameters

u_sm = GM_Mars / (GM_Mars+GM_sun);
u_md = GM_Deimos / (GM_Deimos+GM_Mars);
u_mp = GM_Phobos / (GM_Phobos+GM_Mars);
%[] Mu constant
%    sm --> Sun-Mars system
%    em --> Earth-Moon system%%

DU_sm = 235.34e6; 
DU_mp = (9236.53+9517.58)/2;
DU_md = (23455.5+23470.9)/2;
%[km] Distance from Sun to Mars

TU_sm = sqrt(DU_sm^3/GM_sun);
TU_mp = sqrt(DU_mp^3/GM_Mars);
TU_md = sqrt(DU_md^3/GM_Mars);
%[s] Time unit 
%[] Mu constant
%    se --> Sun-Earth system
%    em --> Earth-Moon system
rm = 3389/DU_sm;
rs = 696340/DU_sm;
figure
hold on
[L1,L2,L3,L4,L5] = librationPoints(u_sm);
Lse = [L1,L2,L3,L4,L5];
for ii = 1:5
    plot3(Lse(1,ii),Lse(2,ii),Lse(3,ii),'bo','MarkerFaceColor','b','MarkerSize',2)
end
[xs,ys,zs] = ellipsoid(-u_sm,0,0,rs,rs,rs,50);
[xm,ym,zm] = ellipsoid(1-u_sm,0,0,rm,rm,rm,50);
Earth = imread('Earth.bmp');
surface(xm,ym,-zm, ...
        'EdgeColor',[255 140 0]/255, ...
        'FaceColor','TextureMap')
surface(xs,ys,zs,...
    'FaceColor',[255 165 0]/255,...
    'EdgeColor','None');
xlim([-1.3,1.3]);
ylim([-1.3,1.3]);
axis equal
set(gca,'FontSize',20) 
end